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PREMIÈRE THÈSE. 



EXPOSE 

DE LA MÉTHODE DE HANSEN 

POim LE CALCUL 

DES PERTURBATIONS SPÉCIALES DES PETITES PLANÈTES. 



APPLICATION. - OOUPABAISON AVEC LA MÉTHODE DES COORDONNÉES RECTANGULAIRES. 



INTRODUCTION. 



Quand on veut déterminer le lieu d'une petite planète dans l'espace, on part 
des éléments elliptiques qui représentent, à une certaine époque prise pour ori- 
gine, la position de l'astre et sa vitesse, et l'on calcule les coordonnées x,y, z, 
résultant de l'emploi de ces éléments; puis les variations &œ, &y, dz de ces coor- 
données, causées par les perturbations. Pour cela, à l'aide des expressions de 
-r^> —rf- et -^— et des formules de quadrature mécanique, on obtient les va- 
leurs numériques de $x, !ty et ds pour les diverses dates considérées. 

Telle était du moins la méthode adoptée jusqu'à ces deux ou trois dernières 
années, et qui semble devoir être bientôt définitivement remplacée par une nou- 
velle, due à Hansen, dont l'exposition figure dans le troisième volume des J/i^ 
moires de la Société scientifique de Saxe, année- 1 SS^. 

Au reste, le Jakrbuch de Berlin, où se trouvent publiés les documents néces- 
saires à ces sortes de recherches, ne communique plus que les données à l'aide 
desquelles s'effectue l'appUcation de la méthode nouvelle. 

■ Ci 
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C.a EXPOSÉ 

Il peut donc y avoir quelque opportunité à exposer sur cette matière les tra- 
vaux de l'astronome de Gotha. Tel est le but du présent Mémoire. 

Par un exemple numérique, nous avons ensuite essayé de montrer que rem- 
ploi des formules de Hànsen simplifie notablement les laborieux calculs d'au- 
trefois. 

Le présent Mémoire est moins une traduction qu'une exposition de la mé- 
thode de Hansen. Nous avons le plus souvent modifié d'une manière simple les 
démonstrations de l'auteur. C'est ainsi, par exemple, que, pour arriver à l'im- 
portante transformation de la page 17, nous avons substitué à Temploi des ima- 
ginaires une méthode géométrique plus claire, et dont l'idée' nous a été suggé- 
rée par une Note de M. Hoplf , de Bonn, insérée au n*" 1353 des Attronomische 
Nachrickten. 

SI. 

Quand on rapporte à un système de plans rectangulaires fixes les coordon- ^ 
nées d'une planète, en tenant compte de la fonction perturbatrice, on obtient 
trois équations différentielles du second ordre en a;, y, z. 

Si l'on néglige la fonction perturbatrice et qu'on appelle Ç, ïj, Ç les coordon- 
nées de l'astre, les équations différentielles deviennent celles du mouvement 

elliptique, et l'on obtient alors Ç", nj, Ç, — • ^ et ^ en fonction du temps et de 

six éléments constants a, b, c, 

La théorie de la variation des constantes arbitraires permet alors de satisfaire 

aux premières équations en prenant pour x, y, z, ~, -£y ^ les mêmes expres- 
sions que pour I, ri, Ç, --» -^» y» à condition d'y considérer les éléments a, b, 

r, ... comme des fonctions du temps. 

Il s'ensuit que, à chaque instant, les coordonnées de la planète et les compo- 
santes de sa vitesse coïncident avec celles d'un point situé sur une .ellipse et 
dont les éléments sont dits éléments osculaleurs de la planète à cet instant. 

Dans notre méthode, nous prendrons pour plans coordonnés un système en 
mouvement, mais soumis à cette condition que les valeurs des coordonnées et 
de leurs dérivées premières, par rapport au temps, aient la même forme dans le 
mouvement troublé que dans le mouvement non troublé. 

Un pareil système a été nommé par IJansen système de coordonnées idéales. 

Le plan mobile des XY, que nous adopterons, sera précisément le plan de.s 
éléments osculateurs de l'astre. 
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DE LA MÉTHODE DE HABSES. C.3 

Des équations nous feront connaître les coordonnées /• et c ( rayon vecteur et 
angle de ce rayon avec l'axe des X), et d'autres les variables qui fixent la posi- 
tion du plan des XY. 

Pour arriver à ces valeurs de r et de v, nous partirons des éléments elliptiques 
à l'origine, et, par un choix particulier de variables, nous obtiendrons très-sim- 
plement nos inconnues. 

Pour conclure la latitude b, la longitude / et, par suite, les valeurs de x, y 
et z rapportées au plan fixe, nous transformerons les formules connues en d'au- 
tres, d'application extrêmement commode, dans lesquelles les grandeurs du 
second ordre seront mises en évidence et pourront, dans la plupart de* cas, être 



Prenons pour.unité de masse celle du Soleil, appelons m la masse de la planète 
troublée, m', m", . . . celles des planètes perturbatrices, k* la constante de l'attrac- 
tion, c'est>A-dire l'attraction de l'unité de masse sur l'unité de masse à l'unité 
de distance; les équations du mouvement de la planète troublée, rapportées à 
des axes rectangulaires fixes passant, par le centre du Soleil, peuvent s'écrire 



i^ + *'(,+„); = *.(, 



(■) 






en posant 

A étant la distance des planètes, m et m', Soit XYZ un autre système de 

coordonnées rectangulaires; on aura, a, «', a" étant les cosinus des angles que 
fait l'axe des X avec les axes primitife, etc., 

i X= «X -y- <Jy -f- e^z^ 

(3) Y=^x+-&y + ^z, 

Nous définirons système de coordonnées idéales tout système tel, que les 
coordonnées de la planète troublée et leurs dérivées premières, par rapport au 
temps, aient la même forme que dans le mouvement non troublé. 

Les cosinus a, a', a', ... étant des (onctions du temps, les dérivées de X, 

Cl. 
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C.4 


EXPOSÉ 


Y, Z deviennent 




dt 


='S-^l-"'ï-^ê-^§'-'^'- 


rfY 
dt 


= Pt-^'l + P'ï--|-rf + .f 


dZ 
dt 


= -'T.^i%*f%^^%-r^;^'% 



et, pour que X, Y, Z soient idéales, il faudra qu'on ait 

Ida da' da' 

dt -^ dt dt ' 



Telles sont les équations générales de condition pour qu'un système de coor- 
données soit idéal. 

Ces équations (4) se réduisent à deux; et, coinine il faut trois conditions 
pour définir un système d'axes, on voit qu'il y aura une infinité de systèmes de 
coordonnées idéales. 

Pour le montrer, je remarque qu'on a 

Ij: = aX + p Y + yZ, 
^ = a'X + p'Y + /Z, 
r = «*X + p"Y-t-7"Z. 

Substituant dans (4}i et posant 

5^ + P'S'+^ï = '-. 



(S) ['i + ^-. 



.!V^„ 



+ !'ï + /^' 



il vient 



; ^Z — /Y = o, 

(7) { '■X — pZ = o, 

' pY — ^X:= O, 
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DE LA MÉTHODE DE HÂSSEN. C.5 

équations dont l'une est une conséquence des deux autres. Or p, q, r sont pro- 
portionnels aux cosinus des angles que fait avec les axes des X, Y, Z l'axe in- 
stantané de rotation du système; donc, en vertu des relations {7), on arrive à 
cette conclusion remarquable : 

Dans an système de coordonnées idéales, l'axe instantané de rotation coïncide 
avec le rayon vecteur de la planète. 

Cette proposition se présente, du reste, immédiatement comme une consé- 
quence des équations (4). En effet, les équations en -^ï -^j -^ expriment alors 
que la vitesse relative coïncide avec la vitesse absolue. La vitesse du point 
coïncident du système est donc nulle à chaque instant, ce qui montre que ce 
point se trouve sur l'axe instantané de rotation. 

La troisième condition étant arbitraire, je prendrai Z = o, c'est-à-dire que le 

plan des XY contiendra à chaque instant le rayon vecteur. On a aussi -^jT = o? 
ce qui prouve que la vitesse relative et, par suite, la vitesse absolue se trou- 
vent également dans ce plan. Ce plan des XY, contenant la tangente à la courbe 
et le rayon vecteur, est donc le plan des éléments osculateurs delà planète 
troublée. 

Par suite de l'hypothèse Z = o, la première des équations (7) donne r = o, 
c'est-à-dire 

lj3rf«-}-p'da'-i-p''d«"=o, 

dont l'une est une conséquence de l'autre. 

On voit alors que les variations des quatre cosinus a, p, a', ^' sont mainte- 
nant liées à celles de a" et JS". 

Combinons, en effet, les deux équations précédentes avec les suivantes : 

«da -h a'dct' -i~ o^d<x"= Oj 



il en résulte, par exemple. 

Or, d'après une propriété bien connue, on a 
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C.6 EXPOSÉ 

d'où 

irf« = 4 t'a", rf& = ^ d&", 
1 7 

7 7 

Il nous faut obercher les équations qui pennettront de déterminer X, Y, a" 

et P". 

I,a fonction iï est exprimée en fonction de x, ^ et z ; cherchons ses dérivées 
jiar rapport à nos nouvelles variables X, Y, Z. 

On a 

.T=aX + pY+7Z, 
,T=<"'X+|3'Y + /Z, • 
j =a'X + ^ï-1-/Z. 

Pour avoir -^i regardons x, ^ et a comme des fonctions de X, et l'on aura 

rfX~*' rfT~^' dz~'^' 

de même ^ = .', ^=P', g = /i 

(ù „ rfa „„ d» n 



et L'on a 
d'où 



da _dtk dx da dy - du di 
dX~ dxdX'^'^ dX~^^ dX' 






et, d'après le système de coordonnées choisi, il faudra ici, après la substitution 
des valeurs de x, y et z, faire Z = o. 
Avec l'hypothèse Z = o, il vient simplement 

x= «X +PY, 
^ = .'X + ?'Y, 
z = «'X + p»Y. 
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DE LA MÉTHODE DE HANSEN. C.7 

Prenons les dérivées par rapport au temps, on aura 

dt di ^*^ dt \ di ^ di ' 

J=:,'!?^+,9'" 0. (,o) X^+Yf .-,-, 

dt dt ' dt ^ ' \ dt dt 

dz .dX ,„dY Iv''*' v*^" 

- ^« --+P —, I X— -(-Y— = 0, 

toujours en vertu de cette proposition que la vitesse absolue coïncide avec la 
vitesse relative. 

Différentions encore par rapport au temps : 





d'x 

~dn 


= " 


rf'X 
dt' 


H-^ 


d'Y 
dt' 


+ 


da. dX 

dt dt 


-S 


dY 
dt 




dt' 


= « 


d'\ 
dt' 


4-P 


d'Y 

dt' 


-t- 


dt 


dX 
dt 


-f 


dY 
Tt 




dH 
dr 


= a 


d'X 

dt' 


+ |3" 


d'Y 
dt' 


-H 


da" 
~dt 


dX 
dt 


-f 


dY 
dt 


'on tirera 


Se. 


d'Y 



















Multiplions la première par «, la deuxième par a', la troisième par «", et 
ajoutons, il viendra 

d'X _ l!f _^ / 1^ H d'z 

<le même 

dt' ^ dt'^^ dt' ^ dt' 

Si l'on multiplie les trois équations respectivement par y, y et y", et que l'on 
ajoute, il viendra 

yun^yiû'^f -dt' = \!dt-^y Tt+^ -d^j-df^xj-di-^'^-t+y in)iû^ 

et, en vertu des équations (9}, remplaçant da., du', dp, d^' par leurs valeurs en 
fonction de da" et rfS", 



/ da' rfX d^ dY\ 

\~d7 Ht '^ Ht lïj 



On a de plu; 
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C.8 EXPOSÉ 

Il est facile maintenant d'obtenir les dérivées des nouvelles variables X, Y, 

a-, |3-. 

Multiplions les équations (i) respectivement par (a, a', a"), (|3, P', ^"), (7, ■/, 7"), 
et ajoutons, il viendra 

/ d'X ,., ,x ,,, , du 
■5?- +*^ (■ + '") r. =*(' + '"' 3x' 

Pour déterminer a" et jS", joignons à la dernière des équations (r i) la troi- 
sième des équations (lo), savoir 

dt dl ' 

de ces deux équations nous tirerons les valeurs de -^ et -^ 



Posant 

il vient 

(.a) 



-S 


-s 


^■%i'^. 


-)i 


-s 


-ï 


»■(.+».) 


-J. 






' dZ 



Les équations différentielles (9), (i i) et (la) déterminent les variables X, Y, 

«,^,7,.... 

A ces variables nous allons en substituer d'autres. 

Au lieu de X et de Y, introduisons les deux coordonnées polaires r et c, r étant 
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DE LA MÉTHODE DE HANSEN. C.9 

le ra^on vecteur de l'astre et v l'angle, dans le plan des XY, de ce rayon vec- 
teur avec l'axe des X positifs, angle compté de o à 360° dans le sens du mou- 



vement. 
On a 



X = r cos i-. 
Y=rsiin' 



-=smi'-r.-{-rcosv- 



lit' dt' dt* dt dl dt' 

d'Y . d'r d'v drdv . rfi* 

—r— =: sint» -;— +rcosi' -7- + 3cos^ -r -7 — ri\a.v -r-- 
dt* dt' dt* dt dt dt' 

Cherchons les dérivées de la fonction perturbatrice par rapport à r et à c, 
on a 

da ^ dû. dH da dY _ da da . 

d? ~ dX'd^ ~^'dY dr ~ dTi'^^^'^ rfr""''' ' 



remplaçant jv^^jf P***" ^^"*^ valeurs tirées des équations (ii), nous arrivons 
aux équations suivantes en r et v : 



(.3) 



Bemarques : 
1" On a 






■*■ "37 — ï ~3r — ^ T. — ~ 



3" Des équations, qui donnent -r- et -rt on tire facilement 
(■4) 



(jX r dv dr 

da 1 da da . 

3v=-t-r T^««''+:i^'"°''■ 
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C.io EXPOSÉ 

Les équations (i 3) détermineront ainsi r elf, que nous substituons à X et à Y. 

Au lie»i de déterminer la position des axes mobiles au moyen des neuf co- 
sinus a, a', ..., nous pouvons choisir d'autres variables en nombre moindre. 

Figurons le plan fixe des œy et le plan des XY à un instant quelconque. 

Soit {^g. i) (j l'angle, dans le plan des XY, de l'axe des X avec la ligne allant 




de l'origine au nœud ascendant; alors v — c sera ce que l'on appelle Vargumenl 
de la latitude. 

Soient b l'angle que le rayon vecteur fait avec le plan des ay; l l'angle que fait 
la projection du rayon vecteur avec l'axe positif des a; (6 et / sont ainsi la lati- 
tude et la longitude de l'astre). 

Enfin,'Soit B l'angle de l'axe des x avec la ligne des nœuds. Alors l — 6 esX 
l'angle de cette ligne des nœuds avec la projection du rayon vecteur. Soit i l'in- 
clinaison du plan de l'orbite sur le plan fixe. 

Le triangle sphérique M £J M' donne 

sîni = sinisiD(c — u), 

. ,, ,, ^ , cosi sinft cos/ i àab 

sinif — 9) ^ taneo -: — = ; -; — r =; — - -; — cosj, 

^ ' " siDf cos6 sini aiib siui 



d'où 
(.5) 
de plus 



r cosi sin(/ — Q) = sm{v — 7) coii, 

j cos£cos(f — e)=co3((' — a), 

[ sini = siiiisiii(i' — a); 

|:r = r cosi cos/, 
y =rcos5siD/, 
z = rsinb. 
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Une fois que l'on aura déterminé r, ç, Q, i et <j, on en déduira /et b, et, pai 
suite, X, y et z. 

Cherchons les équations différentielles en G, i et cr. 

De même qu'il existe une relation différentielle entre les cosinus a, J3; se', ^'; 

«", r. 

- i^a ,, rfa' r a dx" 

il y aura une relation pareille entre $, a et e'. 

Pour la trouver, différenlions les valeurs suivantes de a, a', ce" 

(SE =::cosd cosff + sinff sinacosi, 
a'=: cosffsinS — sinacosScosi, 
[ «"^ — imonai; 



il viendra 



f/f 



dt '^ dt dt ' dt 

rfa' fada M . di 

</( ^ dt ' dt 

Multiplions par ^, ^', P', et ajoutons, il viendra 
d'où 

Telle est l'équation de condition des variables i, 9 et ?. 
Si nous prenons les deux équations 



t sinico89 = ^", 

nous en déduisons, par la différentiation et en tenant compte des équii- 

C.a. 
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lions (13) 
(■7) 



= 3l(Xcosa + Ysinff) — =Xrcoa(i' — a) — i 



- :=Xcot({ — Xsiiio+ Ycosu) -j^ = Xcoti>8in(v — 1) 



Les équations (i3), (16) et (17), prises simiiitanément, déterminent les va- 
riables r, c, 0, i et a. Comptons les constantes arbitraires introduites par l'inté- 
gration. Trois sont données par les équations (16) et (17]. qui sont les valeurs 
de I, $ et 9 au temps f — o, et quatre par les équations (1 3), qui sont les va- 
leurs de r, t>, ^t jt à la même époque. 

On a donc en tout sept constantes, tandis que l'intégration des équa- 
tions f 1) n'en donne que six. Cette circonstance est facile à expliquer, et vient 
de ce que la position de l'axe des X dans le plan des XY est complètement arbi- 
traire; la septième constante déterminant cette position. Comme on peut la 
choisir à volonté, l'hypothèse la plus simple consiste à admettre qu'à l'é- 
poque t = o l'axe positif des X fait avec la ligne des nœuds te même angle que 
l'axe positif de x. Si l'on désigne les valeurs initiales de », Q et<x par („, 80 et ij«, 
la condition énoncée plus haut se traduit par o-^ := $,, ce qui réduit à six te 
nombre des constantes arbitraires. 

Les équations (16), (17) et les deux premières du groupe (i i), dont l'ensemble 
constitue la solution complète du problème, peuvent s'obtenir sur-le-champ, 
d'après la théorie du mouvement relatif, au moyen de la méthode suivante, que 
M. Briot a bien voulu m'indiquer. 

Pour obtenir tes équations différentielles du mouvement relatif d'un corps, 
on sait qu'il faut ajouter aux forces effectives agissant sur le corps deux forces 
fictives, l'une, la force d'inertie du point coïncidant du système, et l'autre, la 
force de Coriolis ou la force centrifuge composée. 

Soient {^g- 3 ) M et M' les positions de la planète aux temps 2 et 2' ; les axes in- 
stantanés de rotation à ces deux époques sont OM, OM'. La vitesse du point coïn- 
cidant M au temps t est nulle; au temps t', elle est la même que si les axes mo- 
biles tournaient autour de la droite OM' avec une vitesse angulaire u'; elle est 
perpendiculaire au plan MOM' et égale au produit de la vitesse angulaire w' par 
la perpendiculaire MP, abaissée du point M sur la droite OM'. La variation de 
la vitesse du point coïncidant pendant l'intervalle de temps t' — t est donc 

m'MP = w'MM' sinMM'P. 
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DE LA MÉTHODE DE HANSEN. C.i3 

Son accélération au temps t est la limite du rapport ta' -y~-- sin MM'P, quand 
l'intervalle de temps tend vers zéro. Elle est égale à wVsinJ, V étant la vitesse 

Fis. 2. 




de la planète; S l'angle de cette vitesse avec le rayon vecteur; elle est perpen- 
diculaire au plan passant par le rayon vecteur OM et la vitesse V, et, par consé- 
quent, parallèle à l'axe OZ, mais en sens contraire. La force d'inertie, pour l'u- 
nité de masse, étant égale et contraire à cette accélération, est dirigée dans le 
sens OZ. 

La seconde force fictive a pour expression a&tVsind*, elle est aussi parallèle 
à OZ, mais en sens contraire. La résultante de ces deux forces est donc paral- 
lèle à OZ et égale à — MVsinî; elle n'aura pas d'influence sur le mouvement 
relatif de la planète dansMe plan mobile XOY, et, par conséquent, les équations 
de ce mouvement relatif seront les mêmes que si le plan était fixe; on obtient 
ainsi tes deux équations 

En projetant toutes les forces sur l'axe OZ, on obtient la troisième équation 



d'où l'on déduit 
L'application des trois premières équations d'Euler, relatives au mouvement 
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Cl 4 EXPOSÉ 

d'un système d'axes rectangulaires mobiles, donoe immédiatemcut: 

— = COti{ — psÎDO-f-^COSff) — (0COttsiD(l ff), 

sini -j- ^ — psiau -\- tj coia = wsin(k — a), 

en désignant par p et q les projections de la vitesse angulaire u sur les axes OX 
et OY ; remarquant que la projection de u sur OZ est nulle et que l'angle a est 
compté ici en sens négatif, à partir de la ligne des nœuds. 
Si l'on remplace u par sa valeur, on a 

tir . -rfn , . 

I rf-î- , .du . , 

I . .(rt . rfû . , . 

\>^ri,~ = Xr — smii'-a). 

Des deux dernières, on déduit 

da .M 

Nous trouvons ainsi rapidement les cinq équations (a) et (c), qui déterminent X, 

Y ou r et V et i, Ô, <3. 

sn. 

Coimaissant c, i, 9 et o-, nous pourrons eo conclure / et 6 à l'aide des équa- 
tions (i5). 

Seulement ces équations ne sont pas les plus simples, parce que les valeurs 
complètes de 6, i et <7 doivent y élre substituées, et la partie constante ne se sé- 
parerait pas de la partie variable introduite par les perturbations. Toutefois, 
une transformation convenable va nous permettre d'arriver à cette séparation. 

Aux formules ci-après, 

/ C0S(C09(/ — ô) :=C08(l' — u), 

1 cosisin (/ — 6) = sîn [v — u) cosi, 
( sîn & = ain I sin (v — a), 
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nous allons démontrer qu'on peut Bubslituer les suivantes : 

/ cos h coa {l — Ot) = cos(f — g) -h AjsinW, 

' cosAsin (/ — 9t) = im{v ^ g) coii^ — AjcosW, 

\ smh = sîa(f — g) sini'o -t-.T, 

dans lesquelles les quantités A et W sont indépendantes de v. 

Pour cela, de l'origine des coordonnées comme centre, décrivons une sphère 
d'un rayon égal à l'unité {^g. 3). Le plan fixe la coupera suivant le cercle x^yo. 




Figurons sur cette sphère la trace du plan de l'orbite osculatrice à l'origine du 
temps, déterminé par $, et l'j, et la trace du plan de l'orbite osculatrice à un in- 
stant quelconque, déterminé par 6 et i. Les deux cercles se coupent en un 
point J. 

Transportons l'axe des ac de Ox^ en OH x. Soit M la position de l'astre ou plu- 
tôt le point où la ligne menée de l'origine à l'astre rencontre la sphère. Mar- 
quons sur le premier cercle un point M«, que j'appellerai ïhomologue de M, et 
obtenu en prenant JMo = JM. Soit de même JE = JX., en sorte qu'on aura 
M«OE = c. Soit JB = JH, et menons le plan bissecteur de l'angle J. 

La droite M«M, étant perpendiculaire à ce plan, a une direction bien dé- 
terminée et indépendante de l'angle v. Appelons W le petit angle HON, cet 
angle étant compté à partir de HO de zéro à 36o degrés dans le sens du mou- 
vement. 

Menons par le centre de la sphère une parallèle 0|x à la droite M»M, c'est- 
à-dire une perpendiculaire au plan bissecteur DNIX. Le plan passant par O/x 
et Os coupera le plan ary suivant une droite OL, située à 90 d^rés de ON, puis- 
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quf ON est perpendiculaire aux deux droites 0/x et Oz. Conséquemment ( 



aura [Jîg. 4)yOL = 'W. Appelons ij l'angle toujours aigu que fait le plan bissec- 
teur avec le plan des xy. 

Si nous joignons N et B par un arc de grand cercle, on aura MB = W, t-t 
l'angle NBG = h. 

Nous pouvons écrire alors 

X = cosbcos(l — 60), 
j = co3&sin(/ — 0„), 
z ^ sin b f 

^0. y,, 2i étant les coordonnées du point M^, on a, en posant HE = g, 

x, = cos(t'— ^), 
j, = sin{i' — g-Jcosi,, 
Zt = sin (c — g) sini't. 

Si nous transportons les axes parallèlement à eux-mêmes {^g. 5) au iioint Mo, 

Fie. 5- 



ï,\". 



on aura, en appelant a;,, y,, z, les nouvelles coordonnées du point M et po- 
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d'où, par substitution, 



DE LA MÉTHODE DE HANSEN. C.17 

M,Q = Z — Z. = *, ung« = A, 
MeMi = Jtaiigv)= A5, 
( Xt= AjsinW, • 
/^, ï= — AjcosW, 



f ï = «, + Z|, 



cosicos(/— e(,) = cos(i' — y} + AssinW, 
cosésiD (/ — On) = sin (v — g}co3i« — AicosW, 
sini = sio (v — g) aia i» -f- s. 

Toutefois, pour arriver aux équattons mêmes de Hansen, qui diffèrent un peu 
de ces dernières, je remarque que, au lieu de considérer un plan fixe dans lequel 
XqH est constant et égal . à 0^, £H ou g étant variable, je puis choisir un plan 
d'inclinaison toujours égale à i,, mais pour lequel £H sera constant et égal à S^, 
tandis que x^H, au contraire, variera. 

Considérons, en effet, le plan mobile XBM, et cherchons à en déduire la posi- 
tion de ce nouveau plan adopté. 

A partir du point X, nous prendrons XB = $a ; nous ferons au point B, avec 
BG, un angle égal à i^; à partir du point N ainsi obtenu nous prendrons 
sur Na7, une longueur NH — NB, et par le point H nous mènerons un grand 
cercle d'inclinaison i,. On voit bien alors que, si l'on prend JE = JX, on aura 
HE = EX = Ô,. 

Il suffit donc, pour avoir les équations résultant de l'emploi de ce nouveau 
plan auxiliaire, de remplacer dans nos considérations précédentes 8^ par G 
(G étant une nouvelle variable) et g par $^. Posant G = Ô, + r, nous obtien- 
drons les formules définitives 

icosA coa (/— e, — r) =co8 (f — flo) -t- Ajsin W, 
cmb sin {/ — 6, — T) = sin(v — ô,) cosi, — A jcob W, 
sîni = sin(v — 6,)aiin"B-t- J- 

Reste à exprimer les quantités A, W et T en fonction de o, $ et i. 

C.3 
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On a 

cosBNG = cos (« — an) = — cosi cosi^ h- sini sini, cos ( a — 9,), 
d'où 

cos an. 
Posons 

K ^ i + cos2)î^ aco5'n= n- cosi cosi, — siDisini, cos ( j — fl»); 

le triangle NBG donno 

Binfsin {» — S,) 

Multipliant par tang i] ou A, il vient 

A sm W = i -^ î 

de même 

cotW8in(ff — Ôt) =co8(ff — fl,)cosi, + siDigCoti. 

Remplaçant sin W par sa valeur, 

C08WsinaM = cosi* sini coa (9 — 9,) H- sini, cosi; 
or 

sinan = aAcos'n; 



d'où 



AcosW: 



G09/|SiRf cos(<r — 9|) + sin/tcasf 



Dans le triangle HBG, on a 

HG = e — 9, — r. 
Posons, pour abréger, 

HG = a et l'angle BHG-C; 
on a alors 

rin» _ Bia(g — fl.) 
sin(/, + C) sinC 

OU 

lina = siii(a — â*) (smi^cotC -l-cosif). 

Remplaçons cotC par sa valeur A cos W, il vient 

sm (, cotC + COSI, = 1 

d'oiî 

(cos/ +cos/,lBÎn(ff — fl,) 
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On a, de plus, 

C0tasill(7 — on) = C08{ff — fl,) COS( 4-8m(C0t(t,-t-C), 

coietsin (o — fl,)sm (», + C) = sm(ij + C) coa(o — 9,)cosi-+-sinico8(j|) -+-C); 

or on a 

sin(i. + C] 



SÎD (iT — 9,j 

d'où 

cosfl 8ÎiiG= 3ia[j(i -(-C)cos (a — 0||)co8i + sinicos(i, + C), 
coset = cosf cos (ff — Sa) (smi»cotC + cosï«) 4- smi(co8i»coiC — smi»). 

La première parenthèse est connue d'après ce qui précède; pour obtenir la 
seconde, remplaçons cotC par sa valeur, et nous aurons 

,, . . cosfir — Stlsîni — rini, 
cosi, cotC — 8ini* = i ■ , 

et, par suite, 

cos(iT— 6,}(i -i-co»(co8/,} — lintsin'a 



Nous arrivons ainsi aux formules suivantes : 

(il)) x= i-+-cosicoBij — sinisini, cos (o — 0, 

8mi£tn(o — 9,] 



AsinW = - 



(ao) 

f Aco8W = 

/ C08fô— ff — D — '^('^~^') (î + coaicos/,) — Mnfainf. 

Examinons les ordres de grande\ir de ces diverses quantités, en considérant 
conune du premier ordre les variations introduites par les perturbations ai, 
AS, Aff. 

Puisqu'on a sini = sini sin (c — a), il en résulte 

j = siiit siQ(c — u) — sintg sin(c — fl,); 
* est du premier ordre. 
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AsinW est de même du premier ordre. 

Ces deux propositions sont, du reste, évidentes géométriquement. 

Pour r, prenons la 'première équation du groupe (i8), et désignons par |3, 
|3', j3", ... des quantités du premier ordre. On peut écrire 

cosi cos (/ — 9, — r) = cos (f — e») -H |3». 

Développons, en négligeant T' et appelant P le coefficient de T, 

Pr = co8(i'— S,:— cosico8(/ — fl,} -t-p», 
Pr = cos { i-— (7 + rr—ûj — cosè cos {/ — e + ô— &J 4- jS' ; 

développant les sinus. et cosinus de {g — Ô,) et (9 — 6,} en séries, 

Pr = cos{v—-3) — cosècos(/— ô) ■+-cos5sin{/— ô) {9—9,) — sm(t' — <T){!T— e,} -H (S'* i 

or on a 

c08(c — a) = cos£cos(/ — 9) 
et 

»m(i' — o)co8i = cos6siu(/ — 9); 

il vient donc 

Pr = sin (t- - a) [(9 — 9,) c08i _ (a _ 9,}] -I- ^\ 
Pr = siii(f-<7)[(6-9,)co3i.-(a-9,)]-+-P"'; 
mais on a 



a — Ba= j dflcosi = {9 — 9,)co8i„ 

(o — 9.) = (9 — 9.) COS.", -H ^, 

cosi, étant une valenr moyenne de cosi entre cosj, et cosi, laquelle ne diflFere 
de cos t'o que d'une quantité du premier ordre ; il en résulte 

Ainsi r est du second ordre 
ïlous allons donner une forme un peu plus commode aux équations (i8V 
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siut sm(a — 60) = p, 
sîiL(Cos((i — lîo) — siai'o = y, 



quandtésdu [vemier ordre. On aura alors 



K = cos(» (cosj + cosig) — qsiait, 
s = qsin [v — So) — pcos {v — 6^); 



AcosW=- 



L(î_j 



d'où l'on déduit 

Aco»W=;tangïo H —-• 

Les équations (18) deviennent alors 

icosAsinf/ — 9, — D =cosi»sin (^— fl,)— i ( Umi-H ? — ), 
' \ «' y^ o xcosf, / 
cosAco8(/ — e,— r) = cofl (t- — M H- 'A 
sini = sini, siii{i' — fl») 4- i. 

Dans le plus grand nombre des cas, les quantités du second ordre F, — > - — 
seront négligeables, et les formules donnant l et b deviendront extrêmement 
simples. 

Nous pouvons dès à présent apprécier l'esprit de la métbode. Pour la solu- 
tion complète du problème, il nous faut connaître, outre r et f , les expressions 
de s, r, p et q. Nous laisserons de côté (, ô et 9, nous bornant à la variable nou- 
velle s. Une fois r, c et J trouvés, aous serons en état, ainsi que nous le démon- 
trerons, de calculer, si cela devient nécessaire, des valeurs suffisammeot appro- 
chées Aep, q et r. 

Voici comment on peut calculer T,^ et y en fonction de r, v et s. 

Reprenons les valeurs de sin et ces ( 9 — S» — T), 



_ (cos/, + cos/) sin [. 



^ = R, 



I _. Ji-|-cosi',cost)cog (g — 9.) — >iat,gini 
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C.aa EXPOSE 

Différentiant, il vient 

coj (8 — «, — r) (<(« — dr) = rfR, 
— !iii(9— «,— r)(<i«_rfr) = <iu, 

d'où 

rfe-<;r = u<iH — Rrfu. 

Or 



cosf \R/ 

Développant la différentielle de t-- et réduisant, nous obtiendrons 

j_ Bia/.sin (»— e,) ,. [sini — sia/,cos(a— S.)] ■ ■. 
dT = ' ai + H sini rfo; 



= coatain(iT — Sft)rfi-H sitiicos(7 — $^)da, 
j = costcos(9 — 6t)di — siQisin(7 — 0») rfo. 



^^^qdp-pd^^ F: 



~J «c< 



la constante étant déterminée de façon que, pour / = o, on ait F ^ o. 

(On vérifie encore ici que T est du deuxième ordre, puisque q etp sont du 
premier ordre.) 

Pour déduire les dérivées de p et g de celles de la fonction perturbatrice, je 
remarque qu'on a 

a." = — sinisinv, 

P''= sinicosiT, 
d'où 

;? = — «"cosfl, — |5"sinfl„ 

y = — 0^ sm6t -i- ^'' cosB^ — sîni». 

Prenons les dérivées par rapport au temps, et remplaçons -jr ^t "j- P*'' leurs 
valeurs tirées du groupe ( i a ) , il vient 

* = Xcosi^r,in(.— «.), 
«'î ^ ■ ^" / û \ 
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Rt les intégrations doivent être prises de manière que p et q soient nuls pour 
/ = o. 

On peut maintenant obtenir une expression simple et commode de F. 

Substituons, en effet, ces valeurs dans l'expression de F, il vient 

'-^^^^ = lrcosig.[y8in(v-0J-pcos(.-9,)], 
cdp— pdq . .da 

et l'on a, avec une exactitude suffisante, 

* ' acos'j, J \''Z/ 

pt\q peuvent aussi être exprimés en fonction de j et de sa dérivée. On a, <'ii 
efifet, 

ysin [y — flj) — pcoi {y — 9,) ^ s-, 

et S étant une coordonnée idéale, puisqu'elle n'est fonction que de h et de v, 
on aura 

De ces deux équations on tire, en remplaçant ^ par -^ [p étant le demi-pa- 
ramètre), 

p = — /co» (f — 9, ) + 7^ T- sin ( i- — fl,). 

Ces équations peuvent être utilisées si les produits ^ et sq donnent quelque 
chose de sensible. On aura approximativement, en posant * = -* 
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§ III. 



Il nous faut maiotenant passer à la détermioation des variables fondamen- 
tales r, V et s. 

Occupons-nous d'abord de r et de v. 

Considérons le plan de l'orbite à l'origine du temps, et soient 

T„, temps'du passage au périhélie; 

pn, demi-paramètre; 

Cot excentricité; 

Rg, longitude du périhélie, 

les éléments de l'ellipse, qui représentent les coordonnées et la vitesse de l'astre 
à l'époque / = o;yétant l'anomalie vraie et s l'anomahe excentrique, on a 
les relations 



I ■*- e,cmf 



tanfi-/=^ 



Si nulle perturbation n'intervenait, ces équations fourniraient les quantités r 
et c à une époque quelconque. Elles déterminent le lieu d'une planète fictive 
circulant dans.une orbite invariable, et dont les coordonnées et la vitesse coïn- 
cident avec celles de la planète réelle à l'origine du temps. 

Par un artifice de calcul, nous allons modifier ces équations de manière 
qu'elles représentent le lieu de la planète troublée dans le plan mobile XY k 
une époque quelconque. 

Pour cela, partant des mêmes éléments, nous remplacerons seulement / par 
une fonction de t, que j'appellerai Ç, et le rayon vecteur qui en résultera, nous 
le multiplierons par (i ■+- v), pour avoir le rayon vecteur réel. La question sera 
ramenée à trouver les dérivées des nouvelles variables Ç et v. 

En mettant Ç au lieu de 2, les grandeurs E,y et r prendront d'autres valeurs, 
que nous désignerons par t], 9 et p. 
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Nous aurons alors, pour le mouvement troublé, 



"^ i-l-c,C0Sf 
; r=p(i + v). 

Beste à déterminer Ç et v. 

Voyons d'abord v. 

De la dernière équation, on tire 

rfc ~ p' {"" do ^ dn ) p> dtY dt '' dt) ' 
Or on a 



et, pour abréger, faisant abstraction du facteur (i + m), 




\ Différentiant de nouveau, 

I p (ù' "f'xdi) ~ 7' dÔ '^ Xe^ T'j \di} ' 

Multiplions la première de ces deux équations par 2 -j-» et la deuxième par r, 
et ajoutons, il vient 

P rf/» "^^ 'dt'^*\f'dt~'' di)~ Fr di ""^7" **/ Wv * 

C.4 
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Substituant dans -j^ nous obtenons 

or les équations (i3} peuvent s'écrire 



■ = *' 
(«) 

I d'où l'on déduit 

l'intégrale étant nulle au temps t = o. 

Substituant ces expressions dans la valeur de -7-^, et réduisant, on a 

* ^' dt' r' p dr r dv r* J dv >" \j '^ J 

Telle est la dérivée de v ; et les deux constantes introduites par l'intégration 
doivent être choisies de façon que, pour f = o, on ait 



Pour trouver la dérivée de Ç, nous avons la relation 

rfy _ *' 

OT d» = âif : donc 

dv k' ,, , dp X'' rfÇ 

et substituant à 3-. sa valeur tirée de la troisième équation {a), 



y Google 



DE LA MÉTHODE DE OAINSEN. 6.27 

Si l'on pose Ç = f + iJÇ, on aura 

Or on a 

^ }Ç sera la perturbation de l'anomalie moyenne, ou, ce qui revient au même, 
de la longitude moyenne. 
Pour trouver la dérivée de la troisième variable s, on a l'équation 

smb = sini(SÎn(f — dg) + s. 

Multiplions les deux membres par r, et posons rs~u, il vient 

u = z — rsint«$ii)[c — 6»). 

Ainsi u est le changement produit par la force perturbatrice sur la coor- 
donnée perpendiculaire au plan fixe. C'est cette variable u dont nous allons- 
chercher la dérivée. 

Diffërentiant deux fois par rapport au temps et introduisant les dérivées de la 
fonction perturbatrice, nous obtenons 

■57r=rfï;-"n'o»>n(*'-9.)^**5:-;îj-"n'.«'«(^-<'«)7Sr- 
Or on a 



1 au M au r u "" H o" 

l — = a" 1- D h 7 ) 

J rfï rfX^*^ d\ ^' d'I 

\ a" = — sinasint, f''=cos<xsini. 
A l'aide de ces trois équations et de celles du groupe {'^), nous déduisons 

d'u i' ,,da u ,, <Al A' da, .. . , . . . , . ,, 

Maison a 

s = sioi sin(v — o) — sîni, aui(n — fl,). 

C.4. 
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C, 28 EXPOSÉ 

Prenant les dérivées et remarquant que f est une coordonnée idéale 

d> < I , u lir dp , . d» ^ dp s'Xàe.nDf dm 

OU, après avou> remplace f par-> — par / (i + v)+p— et ^ par ' - » — ■ 7,» 

on déduira 

imicoa(v — ff)— siQ(,cos(.' — e,) = ^r^— ^p j^u— u^(H-K)pi:ff.aiiiî. 

Substituons dans ~ > nous ol^tenons 

^ ' dt' r' ^ dZ \ f dr r dv j l + 1 d» \dt dH+t/ 

Telle est l'équatioD difTérentielle, d'où nous déduirons la valeur de u. 

Dans les formules précédentes entrent les trois dérivées de la fonction per- 

turbatnce :î-» 3- «t -j=- 

dr dv dï 
Nous allons en chercher les expressions. 
On a, en considérant séparément chaque astre perturbateur, 

da_ m' f, ,\ 



et soit P' l'angle que fait le rayon vecteur de la planète perturbatrice avec le 
plan des XY, on aura 

Z'=r'8in^', 

d'où 

— = "*' A/'sinÛ' 
Pour obtenir j- et -j-t nous prendrons la fonction perturbatrice sous la forme 



il=-^('l — —cosClA, 
l + m\& r" j 

DJgilizedbyGoOt^Ie 
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DE LA MÉTHODE DE HANSEN. C.29 

C étant le lieu de la planète perturbatrice, et P celui de l'astre troublé {fig. 6). 
Appelons XC, v'. 

Fig. 6. 



On a 

ù.* = r*+r'*— arr'^cosCP, cosCP=: cosP'c08((/' — v), 

et, par suite, ' 

rftt m' ( i t \ , ai t I \ m' r 

tir l+m\a.' r"J "^ ^ ' t + mO.' 

Nous obtenons donc les trois dérivées dont nous avons besoin 

= 1 hr'cosB'coslv' — i-) : I» 

(a?) ( -;- = Ar/ cosS' ain f / — c), 

' " j de I + M r \. n 

f <fil m' , , . ,, 

Nous allons maintenant reprendre les trois expressions de ~, -^ et -7^1 et 

leur donner la forme définitive appropriée au calcul. 

Prenons ^r * ^^ différents termes du second membre contiennent tous, soit v, 

soit une dérivée de la fonction perturbatrice qui renferme m' eu facteur. Nous 
pourrons, en ne négligeant que des quantités insensibles, comme nous le ver- 
rons, y remplacer r et p par r„ et f par/, et écrire 

dt'~ /■;""'" 7^ ^ ~ f, dT "*" 7f * J ST ' ' 

et ici r, et /sont le rayon vecteur- et l'anomalie vraie, qui se déduisent des élé- 
ments constants. 



y Google 



noua aurons 



C.3o EXPOSÉ 

Posons, pour abréger, 

dv ' rfr ' rfZ ' 

et remplaçons XI par ~. il viendra, en multipliant les deux membres par le 
carré de l'intervalle de temps adopté w' (pour l'application des formules de 
quadrature), 

<^'" ~ r, rj ^^^ J p,r, /■; 

Nous aurons de même 

d.s\ , rJû j 

dt J dv 

Multiplions par le moyen mouvement fig, et posons 

et, multipliant par u, 

ta — ■ — = Mol — = (1)1 Sa( — a»v (■ 

Quant à u, on a simplement 

— -J. = 2 cos ig J a, 

dou 

*" ■Â^="^^•=°"•-"^;"• 
Tableau des formules en^loyées dans le calcut des pertiahations spéciales. 

R =m'A' I ki-'cosp'coi{v' — c) ^ I» 

S =m'k*hrtr' co9|3'siD{c' — »/), 
Z = m'A'A/-'sinP'. 
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DE LA MÉTHODE DE HANSEN. C.3r 

û)*--=ù)' h2M— — = w JSrfï — — -- w*S r- V, 

-TY = to'^cosi* ~ u. 

"Les formules de quadrature employées sont les suivantes, dues à Encke : 

+ î/.(-J")-à/'(-;")+T^/"('-^)} 

où l'on a 

OU 

Ces formules se trouvent démontrées par M. Airy dans le Nautical Abnanac 
de 1846, et dans l'Ouvrage de Watson. On peut consulter sur le même sujet une 
thèse de M. Gniey. 

M. Puiseux les a également expliquées dans son Cours de la Sorbonne. 

Figurons le plan fixe {Jîg. 7) et les deux orbites dfl la planète et de l'astre 




perturbateur, Jupiter, par exemple ; on a 

XP = f, XE = t^. 
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C.3a 

Or, dans les formules entre (/ 



EXPOSE 
- c, posons 

AP = iv, AE = » 



Il faut calculer »^ et J5'. 

On aura, u' étant l'argument de la latitude de Jupiter, 



DC = u, 


= «'- 


N', 


AD = 


N 


DE = >v' 


-N, 




BD = 


N 


CE = |5'; 










unjtiv' 


-N) = 


= tangl/, cos 


, 


u.6f'= 


tangl 


sm(» 


'-N) 





équations qui détermineront j3' et »/, quand N', N et I seront connus. 

Les données de la question sont (, i", inclinaisons des deux orbites sur le plan 
fixe ; Q, Q' les longitudes de leurs nœuds ascendants. On pourra donc, dans le 
triangle ÂDB, connaissant AB = fl'— fi, les deux angles A = i,B= 180° ~ »', 
calculer N, W et I par les formules 



Il™i(N + N') = .inl(Q'-£3).inlO 

|»ii.ilcosi(N + ri') = ™i(Q'-Q).ml( 

.;Ismi(N-N')=.mj(Q'-a)cosl(, 

!il<»si(N — K')=coil{Q'_n)coii(i 



■■). 
■'+'■), 
■'-■•)• 



Quant à w, on a 

A la distance des deux planètes est fournie par l'équation 

A*^r* + r'* — irr'coi{w' — (v) cos^'. 

Voici donc la marche à suivre : 

Adopter un plan fixe, qui sera, par exemple, le plan de l'écliptique à une cer- 
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DE LA MÉTHODE DE HANSEN. C.33 

taine époque; prendre par rapport à ce plan dans les éphémérides Ci, Q', i et i'; 
calculer alors N, N' et I. 

Les Tables de la planète perturbatrice donnent les longitudes dans l'orbite. On 
en déduira les quantités u', et par suite tv' et J3'. Les mêmes épbémérides 
donnent les valeurs de r'. 

Quant aux valeurs de r, et de /, on les calculera par les formules connues au 
moyen des éléments pris pour base. 

V. 
Application. 

Connaissant les éléments de la planète Eugénie, pour 1857, juillet i, calculer 
l'éphéméride pour décembre 38, et, par suite, les perturbations de quarante eu 
quarante jours dans cet intervalle. 

Voici les éléments adoptés (planètes perturbatrices Jupiter et Saturne) : 

M = 64°.5r.ai'a 1858, janvier o. 

fT = 329.43. 3,t> I Ëquinoxemoyen de1857, 

Q = 148. 5. 3,8 I janvier o. 

i = 6.34.55,8 

(t = i3'io'',955a7 

-1-= 4.43.26.5 

Les formules de quadrature contenant toutes ^ a> pour limite inférieure des 
intégrales, on a pris pour temps zéro la date de juin 1 1 , vingt jours avant juillet 1 . 
On a ainsi obtenu 





Juin II. 


Juillet 31. 


Août 3o. 


Oct. 9. 


Nov. 18. 


Dec. 38. 


Anomalie moy. H... 


ao'.is'.ty'.i 


29! a'.35;5 


iyUd'^hy 


46:37:11^9 


55r24.'3o'i 


64 .'11 .'48 ;3 


■ eic. !.. . . 


«. ..17,5 


3i.3o.33,7 


40. S5. 31,0 


5o. 14.53,0 


59.38.33,7 


68.35.33,3 


. vraie/... 


aî.Sa. 1,68 


34. 4. 9.60 


ji. 7- a.8u 


53.58.43,93 


63.37.47,36 


73. 3.19,14 


— Q+/=.-... 


loS.ag. 2,5 


115.41.10,4 


ia5.44. 3,6 


.35.35.44,7 


145.14.48,3 


.54.40.19,9 




0,400 0788 


0,40» 9530 


0,406 6668 


o,4ii 0738 


0,416 0110 


o,4ai 3i3o 



Éléments des planètes Jupiter et Saturne. 

Jupiter. Saturne. 



: . 1 8 . 40 , 37 
1.55.33,8 



2.39.27,45 
13.28.37,1 
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C.34 EXPOSÉ 

d'où l'on déduit 

Juptler. Satuiue. 

N 189.50.48,70 197. 40 -26,70 

N' 338.56.5a,43 233.12. a,48 

1 5.48.34,64 4.46-5a,7o 

Juio II. JuiUet ai. AoOt 3o. Oct. 9. Nov. iB. Dec. 28. 

w'.... aV-î»'.»»',? a43°ii'.36'i 246'.5o'.a6',5 aSo'aS'.SS'o a54° elss'a 357'.44'.a7' i 

P' +4.a6.ia,5i +4.39.59,48 +4.51.37,10 +5.4.2,4» +5.14.13,91 +5.a3. 6,58 

r'. ... 0,695 569 «.695 95o 0,696 408 0,696 9Î7 0,697 938 0,698 ao4 

Staarne, 

iv' 3i8°5i'.5o',a Sao'ailag's 3ai?5i! 4*4 3a3rao!33'4 3a4?49!56'4 Baôrigîia'g 

p" +4.5.33,90 +4.1.37,05 +3.57-30,58 +3.53.14,70 +3. 48. 49,6a +3.44. i5, 60 

r'. ... 0,955794 0,955978 0,956177 0,956390 0,956618 0,956 859 

Constantes employées. 

m' V = — 7 — E-» m' 5 = cï -pi 

1047,89 35oi,6 

(logtû'A'm'= 4,654967} jr, (log(ù*A*m'=!4,i3ioi6)î), 

log — = a,484o84i logOg = 0,434^70, 

logf/« = 3,898 iSa, loge« = 2,915698. 

\a^\lpl=: o,aiS8o7, 

log~î=;= 1,54861a, 

log -^' r^ 1,922864. 

Type du calcul. 
On n'a inscrit les développements que pour les deus dates extrêmes. 
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/. a3.5a. 1,7 73. 3.19,1 

«/— w 134. 3.26.2 io3. 4. 7,î 

P' +4.a6.ia,5 -i-S.aS. 6,6 

1, sinp' + â,S88 5i4 -i- 2,972 436 

l.r' .-0,695 569 0,698 ao4 

1, r. 0,400 079 0,421 3i3 

I.COS^' -^■î>998 697 -^1,998 »79 

l.cos(»v'— «-) — î,84a aao ~T,354 3Î7 

\.r' 0,800 i58 0,842 6a6 

0,690 129 0,590 980' 0,660 787 0,553 78a' 

i.r" 1,391 i38 1,396408 

l.(r'+r"} 1,490 187 i,5o3 4i3 

0,445 5oi o,5i5i 69a* 0,804 541 o,73o 45o- 

I. arj-'cO8P'C08(w'— «-} i,a37 595 0,77a 963 

1. i» 1,683 096 1,577 5o4 

1.4 0,841 548 0,788 75a 

1. i' a,5a4 644 a, 366 ai6 

l.-î, 3,475 35iS 3,633 744 

o,a4o 843 0,437937' 1,939094 0,271644* 

J. -il 3,913 393 3,905 388 

1. ;■'• a.oSe 707 a, 094 6ia 

l.A -5,716 199 -3,57a 838 

l.r'c08P'coB{*v'— te) —0,536 486 -o,o5o 6ao 

l.Uiig(«^— tv) :. —0,014 a94 — o,634 a67 

lhr'imp'cos(tv'—w) H-â,a5a 685 -1-3,633 458 

0,141 039 0,377 aSo* o,a3o 8a3 o,43i 599* 

1. .^ -^ 3,875 435 -+-2,o55 057 

I.pAr'cpsP'coaCw'-**')-^]- -1-5,016474 -f.3,B54a8i 

I. u'K -t- 6,668 949 — 6,5o6 756 

l.Ar'cosP'aintw'— (c) — 3,a66 979 — â,a57 735 

].hT'cm^iia{i^'—tv] — 2,667 o58 —2,679038 

1. «'S - 5,319 533 - 5,33i 5i3 

m'S — ao8,7o — ai4,54 

1. — -*- 6,368 870 - 6,083 443 



* Nombres serrant an calcul des logarithmea de Gauss. 



C.5. 
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Tablead I. (Suite.) 

18S7. JcpnER. 



Juin r 



r'tànp' :. +T,584 o83 

v'Zcmi, -7,9i9 885 

oi'Zcosi, — 8,91 

~ -6,919 454 

Sio/. 4-T,6o7 044 

vfSdl -:-a,oa6656 

^«/S<A !-i,9495ao 

rj i,aoo a37 

i^'J^ -:o.,49^»3 

iiiîf!^ -F,o,o58. 

/"'■ '■ 

^ 5,475 046 



4-7,670 640 

- 7,893 081 
-7,8" 

— 6,910 aoo 
4- T,98o 714 
-3,008 345 
-a,93i ao9 

i,a63 939 
-7,66? 370 

-7,375 008 

5,411 Hi 



)8»7. Satiikne. 



/, a3.5a. 1,7 

a/—tp 2i3.a2.47,7 

p' +4. 5.33,9 

■ wn^' +3,853 5a3 

.r' 0,95s 794 

,r 0,400 079 

,cosp' -'-■Î1998 891 

.cos(«''-«'} — ï.ff»" 707 

. r' 0,800 l58 

1,143 795 

.r". ï,9ii 588 

,(r'+r") 1,943 953 

0,521 8ia 

. arr'c08P'C08(w' — H') + » ,577 5oi 

. i* -(-a,099 3i3 

.i ",o49 65? 

.i' 3,i48 971 



7Î- 3-19.1 

171.38.53,0 

+ 3.44.15,6 

+ â,8i4 170 

0,956 859 

o,4ai 3i3 

+ ïi999 075 

-1,995 370 

o,84a 636 

1,106 48a 

1,913 718 

1,949 10e 

0,460 a44 

1 ,673 647 

a,i33 891 

1 ,066 946 

3,300 838 



* Nombres servant au calcul des logarithmes de Gauss. 
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Tableau I. (Suite.) 

18S7. Saturnb. . 

Juin 1 1 ■ Décembre 38. 

p ' 5,85i oag S, 799 '^^ 

ï,g6o 307 0,381 S89' o,o5C 618 o,33o a6i * 

4ï 3,i3a 618 3,119 4a3 

r" a,867 38« a, 870 677 

A - S, 81 1 a36 - 5,855 780 

r'c08P'cos(w'— «') —0,876 39a —0,951 3o4 

ttB%[iv'—iv) -t-T,8i8 8o3 — T,i66 766 

Ar'c<»^'C08(«^— w) -^3,687 6a8 +3,807 »M 

o,a38 610 0,436 5ao* 0.456 SgS o,586 609' 

,4 H-5,a5i 108 3,aao 475 

rAr'cOBp'co8(«^- «-) - £,]■ + 3,489 718 + 3,676 870 

«.'R -1-7,600 734 -*- 7,807 886 

Af'C0Bp'Mn(«''-H'} -1-3, 5o6 43i -î,975 84o 

Air'cogP'sintw'— «•) -(- 3,906 5io -3,395 i53 

«•S. -1-6,03? 5a6 ■ — 7,5a6 169 

«'S -1-10,90 -3,36 

^ -i-?,aao655 , -1-7,386 573 

r'tiaP' -1-1,809 3i7 -^-7.77' <»9 

».»ZC08i, —8,748 697 —8,754 953 

«'ZCOS/, • -o,56 -0.57 

ï!? +7,637 447 -7,104 «56 

sin/. ■.-.. -«-î,6o7 044 -i-ï,98o 734 

-JSdt. ..-. ; — o,7"fi 838 •+■ o,aio 586 

l^àufSdl -0,63970a -.-1,133450 

VPt 

r* i,aoo a37 i,a63 gSg 

iîl*^ -ï,439 465 -hî,B69 5ii 

V^ '"• 

<r,gia/*.'S -t-i,7a8 575 -9.569664 

P$ 'V 
^- 3,475 046 â,4ii 344 



* Norobra servant au calcul des lo^rithmes de Ganss. 
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Juin II. . 
Juillet ai. 
Aofit 3o. . 
Oct. 9. . . 
Nov. 18. 
Dec. as. 







Tablea 


u n. 












Calcul de 


^S&dt. 












1857. Jupiter. 










w'S. 










«/Sd(. 


+■ ao9,a8 


—308,70 




+ 5,45 






+ 106,33 


^- o,58 


— iaa,7' 


-14,01 


+ 6,04 


+ 0.59 




- 109,36 


- aai,i3 




- 7,97 




+ 0,45 








-a3o,68 




+ 6,49 




-14 


- 336,o8 


- 451,81 




- i,4B 




-.- o,3i 








-a3a,i6 




-H 6,80 




-i3 


■ - 568, o5 


- 683,97 


-aa6,84 


+ 5,3a 


-1- 6,98 


+ 0,18 




~ 798,13 


- 910,81 


-214,54 


-i-ia,3o 


+ 7,«3 


-+- o,o5 




—1019,40 


-iia5,35 




+19,33 











Juin II 










-t- 0,09 




- a, 19 


Juillet ai . . . 




+ 8,71 






-t- 8,80 




-a,69 


Aoflt 3o 




-.- 6,0a 






-H 14, 8a 




- 3,01 






+ 3,01 






+ 17.83 




- 3,18 












■+■ 17,66 




-3,19 


Dec. a8 




- 3,36 






-H i4,3o 




- 3,o5 



— 0,67 

- o,5o 



5,3i 

■ 4,65 

■ ia,a9 
' 16 ,58 



yGqoglc 




C.39 





- '.n 


o,te 


- ',7) 


•,9< 


~ 10.35 


I,i3 


- >9.39 


0,3a 


- 5j."3 




-91,85 
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TABLEiu II. (Suite.) 
Calcul de m. 



JUÏD II... 

Juillet ai.. 
Août 3o. . . 
Oct. 9.... 
Nov. le. . 
Dec. a8... 



- o.m 




- 0,410 








+ o,-m 




— 0,175 


- o'o7 


H- 0,034 




- 0,585 








— 0,319 




- 0,076 


- 0,08 


-o,r.sS 




- o,6«i 








-0.980 




-t- 0,00a 


-0,73 


- i,«5 




- o,659 








- 1,659 




+ o,i35 


-a,o5 


- »,9I4 




-o,5K 








- a,i63 




+ o,a57 


-3,69 


- 5.077 




- o,.67 








-;.,43o 


» 


- o,Î77 


-6,»9 


- 7,So7 




-4- 0,110 







Clierchons la position de la planète pour la date 1857, décembre a8. 
On a (unités du septième ordre) 



On trouve alors 



v= — 59,38, « = —91,85, ÎM=— 6*,a9. 

T = 73»3'ia^78, 

(f -t- it, = i' = 3o3"45'i6',38. 



Connaissant 9 et a, on aura les valeurs de / et de 6 par les formules 

/ cosbcQs{l — 9,) =cos(i' — S»), 

1 GOE&sm(/ — Ôj) :=sin(i' — 6a)cosit — tuig>t -• 



f siuh 

l =3o2"54'*' 



= sin(*' — Ofl) sini'o 



Stmb = + a, 690535a, 
cos4 = + 7,999477*> 
= p(H-v) = (log) 0,43 i 3094; 
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d'où (coordonnées héliâcentriques écliptiques) 



:c=:;-cosicos/, 


/ log« + 0,1557345, 


j = rcosi9in/, 


logj-- 0,344 8656, 


.=,si„i; 


( legs + 7,111 8346. 



Passant de là aux coordonnées héliocentriques équatorialei (e = tnriinaison 
<le l'écliptique}, 

X| = x, f logjr, -f- 0,1557345, 

/, =jco9ï — z sins, I iogy, — o,3i83843, 

Xi =^sia> -»-2C0Si; ( log^i — 1,8819733. 

Par la méthode des coordonnées rectangulaires, on a trçuvé 

Dec. a8. 3j:, = — î^i, *;•,=— a43, iz,^ — laa; 

d'où 

/ If^o-, = + 0, 1557337, 
î log^i = — o,3i8a85a, 
( loga, =— 1,8819734- 

L'accord est parfait. 

Les équations différentielles qui servent à établir les variations des coor- 
données dans la deuxième méthode sont les suivantes : 

-Sr = '«' (-^ - 7:) + 7j (/? — ^=h 
dans lesquelles on a 
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f=H 



-5j, 



..3.4 ■' 



L'examen des calculs à effectuer dans les deux méthodes est tout k l'avantage 
de la première. 

Il n'v a de différentielle un peu longue à calculer que celle de v. Celles de u 
et de Î5I sont extrêmement simples et s'obtiennent avec rapidité; tandis que, 
dans le second procédé, il y a trois intégrales toutes pareilles à effectuer. D'aulre 
part, les variatious âx, ôy, àz augmentent plus rapidement de beaucoup que 
les nouvelles variations adoptées, ainsi qu'on peut en juger par les résidtats com- 
paratifs suivants ; 

/ V =-59 iS.r^~H', 

< u ~ — 92 ■ ^ = — ^4^1 

((ÎM = — 6,3; (53=— I2Î. 

11 nous a été ainsi permis d'agrandir l'intervalle de temps et de le fner à 
quarante jours, tandis qu'ordinairement il n'est que de trente. 

Dans l'expression de j4< par exemple, nous avons mis, àla place de p, rets, 
partout 7*0 ety, ce qui est absolument permis pour un temps assez long; mais, 
une fois le calcul en marche, rien n'empêche de recourir aux véritables valeurs. 
En effet; on connaîtra alors un certain nombre de valeurs de v et de 9fll, et 
rextraj>olation fournira des valeurs de v et de âM suffisamment précises pour 
le calcul ultérieur de p, ret 9. 

On a refait le calcul pour la dernière date avec les valeurs rigoureuses de p, 
/•et f, sans trouver de modification au calcul précédent. Ajoutons que dans notre 
travail nous nous sommes servis avec beaucoup d'avantages des Tables de loga- 
rithmes de fiauss, d'addition et de soustraction, et de l'édition à sept décimales 
de Wittstein, 

Il y a, dans l'emploi de ces Tables, que l'on ne saurait trop, suivant nous, re- 
commander, une abréviation notable du travail et un moyen d'éviter les erreurs. 

Il nous reste à dire quelques mots des termes qui ont été négligés dans les 
expressions des dérivées des variables, à montrer que ces termes sont extrême- 
ment petits, et à faire voir en outre que, dans le cas où ils atteindraient tme 
valeur sensible, on pourrait facilement les obtenir au moyen des tableau^ de 
calculs déjà écrits. 
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DE LA MÉTHODE DE HANSRN. C.43 

Dans 04* -r^ nous avons omis le terme 

on peut l'écrire 

or w/SAse trouve au tableau (II), et/),/* au Ubleau (I). 

En prenant pour unité la septième décimale, dans notre exemple, le dé- 
cembre 28, ce terme n'atteint pas 0,01. 

a** Dans u* -j-^i nous avons négligé d'abord la partie 



que 1 on peut écrire 



I^a parenthèse est connue d'après le tableau (I) et u d'après le tableau (11). 
Ce terme est d'une extrême petitesse. 



, * da du _ w' „ 

u*5 estdans le tableau (I); <ù y- se calcule au moyen des nombres du tableau (II). 

C:e terme est encore absolument négligeable. 

Quant au dernier X 3- « -7" » il est du troisième ordre. 

^ dv dt i-hi 



Les parties négligées Sfi^ «v*, — a -^j= w / S(7( x v sont complètement in- 
sensibles; car les quantités 3^^ uv et — ^) «> / ^dl (dans notre exemple), dont 
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(:.ij4 EXPOSÉ DE LA MÉTHODE DE HANSEN. 

la pl«s grande n'atteint pas 5 secondes, sont multipliées par v, qui est plus petit 

que 0,000006. 

Par des considérafions toutes pareilles, on arriverait à apprécier et à évaluer 
au besoin les variations de 5 et de i. 

Quant aux expressions sp et sg, ainsi que F, on voit, par les formules des 
pages 35 et 36, que leurs valeurs sont insensibles. 



Vu el approuvé. 
Paris, k 1 1 jaovier 1877. 
Le Dotem de là Faculté des Scierces, 
MILNE EDWARDS. 
Permit d'imprimer. 

l.E Recteur de l'Académie de Paris, 
A. MOURIER. 
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